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ABSTRACT 

We study a complete, colour-selected sample of double-degenerate binary systems 
containing extremely low mass (ELM) ^0.25 Mq white dwarfs (WDs). We show, for 
the first time, that Milky Way disk ELM WDs have a merger rate of approximately 
4 x 10~ 5 yr _1 due to gravitational wave radiation. The merger end-product depends on 
the mass ratio of the binary. The ELM WD systems that undergo stable mass transfer 
can account for >3% of AM CVn stars. More importantly, the ELM WD systems that 
may detonate merge at a rate comparable to the estimated rate of underluminous SNe, 
rare explosions estimated to produce only ^0.2 Mq worth of ejecta. At least 25% of 
our ELM WD sample belong to the old thick disk and halo components of the Milky 
Way. Thus, if merging ELM WD systems are the progenitors of underluminous SNe, 
transient surveys must find them in both elliptical and spiral galaxies. 
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1 INTRODUCTION 

We use the first complete, well-de fined sample of ex tremely 
low mass (ELM) «;0.25 M WDs (|Brown et alJl201Ch to cal- 
culate the space density and merger rate of these unusual 
double-degenerate systems. ELM WDs are noteworthy be- 
cause 1) they comprise <0.2% of a ll spectroscopically con- 
firmed WDs (jEisenstein et aJj|2006| ) and 2) the Universe is 
not old enough to produce ELM WDs through single star 
evolution. Radial velocity follow- up of the 12 ELM WDs 
in th e Hypervelocity Star Survey ()Brown et al.ll2005l . [2007I . 
2009) reveals that they are all binaries with short (1-14 
hr) orbital periods. Remarkably, six of our 12 ELM WDs 
will merge in less than a Hubble time. A simil ar merging 
WD system (J1257+5428, iBadenes et ail |2009h generated 
excit ement when its compan ion was thought to be a neutron 
star (|Thompson et ai1l201Qj ). Improved s pectroscopy shows 
that the system contains a pair of WDs IjMarsh et alj|2010l : 
iKulkarni fc van Kerkwiikll201Ch . just like the ELM WD sys- 
tems studied here. 



* Based on observations obtained at the MMT Observatory, a 

joint facility of the Smithsonian Institution and the University of 

Arizona. 

f Spitzer Fellow 



Depending on the stability of mass transfer in the Roche 
lobe overflow phase, ELM WD binary mergers produce two 
possible outcomes. ELM WD binaries with mass ratios <0.2 
experience stable mass transfer and form AM Canum Ve- 
naticorum (AM CVn) systems, a class of ultracompact bi- 
naries i_thjitconsist of a WD accretor and a helium donor 
star (|Warned 1 19951 : iNelemans! 120051 : ISolheiml |2010j ). With 
their short <1 hr orbital periods, A M CVn stars are im- 
portant gravitational wave sources l|Hils fc Bender! |200d : 
INelemans et al.ll2004l ). Three proposed formation channels 
for AM CVn stars involve three different types of donor 
stars: WDs , helium stars, or slightly e volved main-sequence 
stars (e.g. IPodsiadlowski et al.l 120031 ). Distinguishing the 
dominant formation channel for AM CVn stars is difficult, 
however, because all three channels lead to the same helium 
mass-transfer state. Studying the space density and merger 
rates of progenitors such as our sample of ELM WD binaries 
places an important constraint on AM CVn formation. 

ELM WD binaries with mass ratios >0.7 experience un- 
stable mass transfer that will result in a merger or possibly 
an explosion. Whether ELM WD binaries merge or explode, 
and how they do so, are open theoretical problems. This 
Letter is focused on observations: we compare the observed 
merger rate of ELM WD binaries with one possible outcome 
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relevant to current transient surveys, underluminous super- 
novae (SNe). 

Underluminous SNe, such as SN 2005E and SN 2008ha, 
are rare types of supernova explosions that are 10-100 times 
less luminous than a normal SN Type la and have only 
~0.2 Mq worth of ejecta. One proposed origin for under- 
luminous SNe is the detonat ion of a helium layer on a 
sub-Chandrasekhar mass WD (|Nomotolll982l ; IWooslev" et al.l 
|1986T >. iGuillochon et al.1 (\20m argue that instabilities in 
the accretion stream can detonate the helium layer. In the 
".la" scenario, a sufficient mass of helium accumulates in 
the AM CVn mass transfer phase to ignite and explode 
l|Bildsten et al.ll2007l ; IShen et al.ll2010T ). It is unclear which 
scenario actually occurs for a merging ELM WD, however, 
because the physics of the merger process, such as the mass 
transfer rates, are poorly constrained. 

Observationally, ELM WD merger models can suc- 
cessfully explain the stellar environment, nucleosynthesis 
products, spectra, and light curves of observed underlu- 
minous SNe. Half of the known calcium-rich underlumi- 
nous t ype Ib/c SNe are o bserved in old-population environ- 
ments l|Perets et al.ll2010T ): SN 2005E is located 11 kpc (pro- 
jected) above the disk of its host galaxy, while SN 2000ds, 
SN 2005cz, and SN 2007ke are located in elliptical galaxies. 
These objects cannot be explained by the core-collapse of 
massive stars, but are easily explained by merging WDs. A 
low ejecta mass of calcium-rich material is a n atural result of 
burn ing helium-rich material on a WD (e.g. IWooslev et al.l 
Il986h . The rapid light curve decay of SN 2010X, for ex- 
ampl e, can be explaine d with the .la model l|Kasliwal et al.l 
|2010| ). IWaldman et al.l i|2010l ) are able to model both the 
observed spectrum and the light curve of SN 2005E with 
the detonation of a 0.2 Mq He-layer on a 0.45 Mq WD. 
Half of our merging ELM WD systems have similar mass 
ratios. Measuring the rate of our ELM WD mergers thus 
provides an additional observational constraint on links to 
underluminous SNe. 

Our Letter is organised as follows. In Section 2 we dis- 
cuss the ELM WD sample. In Section 3 we derive the space 
density of ELM WDs. In Section 4 we estimate the merger 
rate of ELM WDs in the Milky Way, 4 x 10~ 5 yr -1 . This 
first estimate suffers from factors of 2 uncertainties, however 
it demonstrates that merging ELM WD binaries are viable 
progenitors of both AM CVn stars and underluminous SNe. 



2 SAMPLE 

Our sample of twelve ^0.25 Mq ELM WDs comes from a 
radial velocity survey of stars selected uniformly by mag- 
nitude 15 < g < 20.5 and colour. The approximate colour 
limits are 0.1 < {u - g) < 0.9, -0.45 < (q - r) < -0.2, 
and -0.5 < (r - z) < 0; iBrown et all (|2007l . 120091 ) provide 
the exact prescription. The colour-selection was designed to 
find halo stars with late-B type colours. However, 15% of the 
survey stars are WDs, 12 of which have ^0.25 Mq. Follow- 
up spectroscopy reveals that eleven of t he 12 ELM WDs ar e 
binaries with 1-14 hr orbital periods (Bro wn et al.l 12010). 
The single non-variable ELM WD is consistent with the 
number of pole-on systems expected in a sample of 12 non- 
kinematically selected targets. Thus we assume that all of 
the ELM WDs are binaries, and that the orbital parameters 



of the 11 ELM WDs fairly sample the underlying distribu- 
tion. 

We estimate luminos ities and dis t ances for our ELM 
WDs using the updated iPanei et al.l (120071 ) evolutionary 
tracks for He-core WDs. The luminosity estimates encom- 
pass +8.0 < M g < +9.9, and co rrespond to helioc entric 
distances of 0.3 kpc < d < 3.1 kpc (|Brown et al.ll2010l ). The 
sample covers a 9800 deg 2 region of the SDSS Data Release 
7 imaging footprint, and is currently 90% complete in this 
area. Thus we assume that there are 10% more ELM WDs 
in our colour-magnitude range when estimating the ELM 
WD space density. 



3 DENSITY 

The simplest approach to estimating the space density of 
ELM WDs is to divide the observed number of ELM WDs 
by the survey volume. This simple approach provides a ro- 
bust lower limit, but only for ELM WDs sampled by our 
survey. ELM WDs are missing from our sample 1) because 
they cool and 2) because they merge. Thus we need to cor- 
rect our survey for ELM WDs that have cooled outside of 
our colour range and that have merged. Our survey samples 
ELM WDs that formed in the last 1 Gyr, thus that is what 
we consider here. For a constant formation rate, we will show 
that older ELM WDs contribute only a few percent to the 
present merger rate. Because the contribution of these de- 
layed mergers is much smaller than our uncertainties, the 
unknown formation history of ELM WDs does not signifi- 
cantly affect our conclusions. In this section we estimate the 
present number of ELM WDs that exist in the Galaxy. In 
the next section we estimate the ELM WD merger rate. 



3.1 Cooling Time Correction 

Once formed, WDs cool and become redder and fainter 
with time. The rate of cool i ng dep ends sensitively on WD 
mass. Updated IPanei et al.l |2007) evolutionary tracks for 
He-core WDs show that 0.17 Mq WDs spend ~1 Gyr cool- 
ing through our colour (effective temperature) range with 
approximately constant luminosity, M g ~ +8.0. On the 
other hand, 0.2 - 0.25 M WDs spend only ~0.25 Gyr cool- 
ing through our colour range and with declining luminosi- 
ties, M g dropping from +8 to +10. The difference between 
the 0.17 M Q and 0.2-0.25 Mq He-core WDs is explained 
by the thermo nuclear flash thres hold in the hydrogen shell 
burning phase (|Panei et al.l 120071 ) . We note that the tracks 
have the WDs spending 10 7 — 10 8 yr at hotter tempera- 
tures, or 1%-10% o f the time spent in our colour range. 
ISerenelli et alj (|200ll ) tracks yield similar numbers. 

Clearly, our colour selection samples ELM WDs that 
have formed within the last 1 Gyr. In addition, the relation 
between mass and cooling time skews the number of ELM 
WDs of different mass that we observe. Each observed WD 
implies the formation of (1 Gyr)/t t s such objects in the 
last Gyr, where t & s is the time spent cooling through our 
colour range in Gyr. We take t b s = 1 Gyr for the eight 
~0.17 Mq WDs in our sa mple, and t i, s = 0.25 Gyr for the 
four 0.2 - 0.25 M WDs j|Serenelli et all 1200 ll ; Fanei et all 
l2007h . This correction factor links the observed sample of 
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Table 1. ELM WD Properties 



ID 



Mi 
(M ) 



^obs 

(Gyr) 



Mi/M* 



(Gyr) 



J0755+4906 





17 


1 





151 


0.17 


5.8 


J0818+3536 





17 


1 





513 


7.31 


1.0 


J0917+4638 





17 


1 





489 


30.3 


1.0 


J0923+3028 





23 


0.25 





528 


0.11 


9.4 


J1053+5200 





20 


0.25 





614 


0.13 


7.6 


J1233+1602 





17 


1 





142 


1.69 


1.0 


J1422+4352 





17 


1 





309 


34.7 


1.0 


J1439+1002 





18 


0.25 





292 


44.1 


1.0 


J1448+1342 





25 


0.25 








1.0 


J1512+2615 





20 


0.25 





555 


140.3 


1.0 


J1630+2712 





17 


1 





242 


12.5 


1.0 


J2119-0018 





17 


1 





163 


0.42 


2.4 



60° 



ELM WDs to the total number that have formed in the last 
Gyr that presently exist in our survey volume. 

The relation between luminosity and temperature also 
affects the volume over which we observe different mass ELM 
WDs. This is accounted for in our 1/Knax density estimate 
below, with the implicit assumption that the population of 
ELM WDs is the same in all volumes. 



3.2 Merger Time Correction 

ELM WDs are missing from our sample because they merge. 
The merger time due to angular momentum loss from grav- 
itational wave radiation is 



(Mi + M 2 ) 
MiM 2 



1/3 



P 8/3 x 10 7 yr 



where the masses are in Mq and the period P is in hours 
l|Landau fc Lifshitzlll95sh . 

Both the mass ratio M1/M2 and the merger time r de- 
pend on the mass of the unseen companion and the unknown 
orbital inclination (M2 sini). An edge-on orbit with i = 90° 
gives the minimum possible companion mass and the longest 
possible merger time. The average inclination for a random 
stellar sample is i = 60°. We assume that the ELM WDs 
have the average inclination in the following discussion and 
in the values tabulated in Table 1. 

Four of the ELM WDs have merger times less than the 
time spent cooling through our colour range. If the forma- 
tion rate is constant over the last Gyr, then the presence of 
these short merger time objects implies that (1 Gyr)/r such 
objects formed over the last Gyr. This correction factor links 
the ELM WDs that presently exist to the total number that 
have formed in our survey volume in the last Gyr. 



3.3 Local Space Density 

We now solve for the total number of ELM WDs that formed 
throughout Galaxy in the last Gyr. First, we combine the 
correction factors such that we do not over-correct systems 
that merge faster than they cool. 



Next, we calculate the local sp ace density of ELM WDs 
us ing the modified 1 /Vin ax method (|Schmidtll975l ') described 
bv lRoelofs et al l (|2007aT ). This method accounts for the vari- 
ation of space density with dista nce above the Galactic plane 
using a disk model. We use the lRoelofs et aL (l2007al) thin- 
and t hick-disk model combined with the iNelemans et al.l 
exponential radial scale length 



P(Z,R) 
Po 



0.98 sech(^) +0.02 sech{-^) 



exp( 



8.5 - R 
2.5 ' 



N cor — max(- 



-)• 



where Z is the vertical distance above the Galactic plane in 
kpc and R is the radial distance along the Galactic plane 
in kpc. Although the thick disk fract ion adopted here is low 
compared to more recent work (e.g. lJuric et al.ll2008h . this 
facilitates a fair comparison wit h the AM CVn population 
derived by dRoelofs et alj|2007al ). 

To solve for the local space density of the disk, we ex- 
clude halo ELM WDs from our sample. Two probable halo 
objects are J0818+3536 and J1422+4352, which have sys- 
temic radi al velocities of —20 2 km s _1 and —193 km s _1 , re- 
spectively. |KiiiceEaL| (|2010r i argue that J1053+5200 is also 
a halo object based on its proper motion. Reliable proper 
motions are unavailable for the majority of our ELM WDs, 
however. We choose to make a self-consistent cut on systemic 
radial velocity. The inclusion of a spurious halo object in our 
disk sample would cause us to over-estimate the local space 
density by ~10%. 

The volume sampled by our survey varies with both ra- 
dial distance R and vertical distance Z. We compute V m ax 
by starting with the SDSS stripe definitions, converting 
them to Galactic coordinates, and then integrating the vol- 
ume over our 15 < g < 20.5 apparent magnitude range for 
each WD absolute magnitude. 

Given the assumed disk model, the completeness- 
corrected local space density of disk ELM WDs is po = 40 
kpc" 3 . The error is dominated by the factor of 2 uncer- 
tainty in our correction factors. Integrating the disk model 
over < R < 25 kpc and < \Z\ < 10 kpc, we estimate 
that 5 x 10 4 ELM WDs formed in the Milky Way disk in 
the last Gyr. 



4 MERGER RATES 

The merger times are encoded in the orbital period distri- 
bution of the ELM WDs. In principle, we can calculate the 
initial period distribution required to produce the observed 
period distribution under the assumption of gravitational 
wave radiation energy loss. Given our small number statis- 
tics, this approach is fraught with uncertainty. 

The approach we take is to consider the observed sam- 
ple of ELM WDs a snapshot in time of the whole population. 
We use the observed orbital periods as a proxy for the under- 
lying period distribution. For reference, Figure [T] plots the 
mass ratio M1/M2 and merger time r for our sample of ELM 
WDs, excluding the probable pole-on system. Filled squares 
in Figure Q] are the most probable values for i — 60° , and 
lines show the possible values over the 90-percentile range 
25° < i < 90°. 

Four of the ELM WD systems in our sample have 
merger times less than 1 Gyr. The short merger times of 
these systems suggest that they represent 70% of the ELM 
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0.75 



0.6 



0.45 



0.3 



0.15 




Figure 1. Gravitational wave radiation merger time versus mass 
ratio for the sample of ELM WDs, excluding the probable pole- 
on system. Mi is the observed ELM WD and M2 is the unseen 
binary companion. Filled squares show the most probable values 
for i = 60°, and lines indicate the possible range of values over 
the 90-percentile range 25° < i < 90°. ELM WD binaries with 
mass ratios below the dashed line will experience stable mass 
transfer and become AM CVn stars; binaries above the dashed 
line will experience unstable mass transfer and may explode as 
undcrluminous SNe. 



of the Milky Way would also increase the merger rate esti- 
mate: star count models suggest th at the stellar halo con- 
tains 10% of the stellar dis k mass (lJuric et alj|2008h . and 
the bulge a similar amount (jWidrow fc Dub inski 2005j). All 
of these corrections act to increase the merger rate, thus our 
estimate is likely a lower limit to the true ELM WD merger 
rate. 



4.1 Comparison with AM CVn Rates 

The ELM WD systems with mass ratios <0.25 will likely 
become stable mass-transfer AM CVn systems. Based on the 
mass functions, the companions to these objects must have 
>0.5 Mq and thus are probably normal C/O core WDs. 
Although the stable mass-transfer ELM WDs comprise 50% 
of the observed merger systems (2 of 4 with r < 1 Gyr, or 3 
of 6 with t < 8 Gyr), the correction factors suggest that they 
account for 33% of the overall ELM WD merger rate. The 
different percentages illustrate the small number statistics 
involved. Thus we consider the range of rates, (1 — 2) x 10~ 5 
yr -1 , with the reminder that our merger rates are likely 

und erestimates . 

iRoelofs et all |2007al lbl) make an order-of-magnitude es- 
timate of the disk AM CVn birth rate of 5 x 10~ 4 yr" 1 based 
on samples containing 5 and 6 AM CVn stars. Comparing 
with the stable mass-transfer ELM WDs in our sample sug- 
gests that ELM WDs contribute at least 2%-4% of the AM 
CVn population. Although this is an order-of-magnitude 
comparison, our observations show for the first time that 
merging ELM WD systems are a likely source of AM CVn 
stars. 



WD systems that formed in the last Gyr, as determined 
by the correction factors calculated above. Combining the 
above numbers, the merger rate of ELM WDs in the Galaxy 
is approximately 0.7 • 5 x 10 4 /1 Gyr ~ 4 x 10~ 5 yr" 1 . We 
get a similar number if we divide the number of ELM WDs 
by the harmonic mean of the observed merger times. 

A 4 x 10~ 5 yr -1 merger rate is likely an underestimate 
for a few reasons. First, some ELM WD systems formed 
long ago and are now merging. In our corrected sample, 5% 
of ELM WDs have mergers times 1-8 Gyr. If the formation 
rate of ELM WD systems was c onstant over the 8 Gyr age 
of the disk (jLeggett et al ] |1998T ). then the contribution of 
delayed mergers is approximately 10 -6 yr -1 . The contribu- 
tion increases if we account for the increased star formation 
rate at early time, but deceases if we account for the delay 
to form ELM WD systems. Fortunately, the contribution of 
delayed mergers is so small that knowing the detailed for- 
mation history of ELM WDs is unnecessary; the ELM WD 
merger rate is dominated by short merger times. 

Secondly, our observations may be missing ELM WD 
systems with merger times < 10 s yr. For a 1 Gyr lifetime in 
our survey region, we expect ELM WDs with merger times 

< 10 8 yr to contribute < 10% to our sample. Objects with 

< 10 s yr merger times have < 1 hr periods, periods in the 
realm of observed AM CVn systems. The contribution of 
ELM WD systems formed with < 1 hr periods may be better 
estimated from AM CVn rates. A single 10 7 yr merger would 
double our estimated merger rate. 

Finally, accounting for the halo and bulge components 



4.2 Comparison with Underluminous SNe Rates 

The ELM WD systems with mass ratios of ~0.55 will 
likely become unstable mass transfer systems. These sys- 
tems include the two shortest merger time objects; the 
correction factors imply that unstable mass-transfer ELM 
WDs account for 67% of the ELM WD merger rate. Given 
the observed mass ratios, there is a 40% probability that 
the ELM WD companions are C/O WDs (M^0.45 M ). 
That means 40% of the unstable mass transfer systems 
may lead to surface H e detonation on a C/O WD (e.g. 
iGuillochon et alj|20ich . Stable- mass transfer system s may 
also detonate as ".la" SNe (e.g. iBildsten et al.ll2007l ). Thus 
the rate of ELM WD mergers that may become underlu- 
minous SNe is (2 — 3) x 10~ 5 yr -1 , with the reminder that 
these rates are likely underestimates. 

For comparis on, the Type l a SNe rate of the Milky Way 
is 5 x 10" 3 yr" 1 (|Li et alj|201oh . ELM WDs are an unlikely 
source of Type la SNe because of their low system mass. 
The relevant comparis on is with underluminous SNe. 

iFolev et al] <|2009t > estimate that a SN 2008ha-type ex- 
plosion is observable to a distance of 40 Mpc, a volume that 
has had 60 known S Ne Type la in the p ast 10 years. Thus, 
based on one object, IFolev et all (|2009h estimate that the 
underluminous SNe rate is ^2% of the Type la rate, or 
~1 x 10 -4 yr -1 . This is a few times larger than the rate 
of ELM WD mergers that may become underluminous SNe. 

We close by noting that at least two of our ELM WDs 
are halo objects and an equal or greater number are proba- 
bly thick disk objects. The halo and thick disk are old stellar 
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populations. If ELM WDs are the progenitors of underlumi- 
nous SNe, underluminous SNe must be found in both ellip- 
tical and spiral galaxies. 



5 CONCLUSIONS 

We study the space density and merger rate of double- 
degenerate binary systems containing ^0.25 Mq ELM WDs. 
This work is motivated by the first complete, well-defined 
sample of 12 ELM W Ds fortuitously targ e ted by the H yper- 
velocity Star Survey (|Brown et all 12003 . 12001 120091 ) . The 
ELM WDs are all cons istent with being short-period bina- 
ries (|Brown et alj|20ich . Remarkably, 6 of the 11 systems 
with good period determinations will merge in less than a 
Hubble time due to gravitational wave radiation. 

The merger rate of disk ELM WD binaries in the Milky 
Way is approximately 4 x 10 -5 yr _ . Although the uncer- 
tainty in this estimate is at least a factor of a few, the 
observations establish that merging ELM WD binaries are 
viable progenitors for several interesting phenomena. ELM 
WD systems that undergo stable mass transfer contribute at 
least 2%-4% of the observed AM CVn population. More im- 
portantly, ELM WD systems that may detonate merge at a 
rate comparable to the observed rate of underluminous SNe. 
Other possible outcomes include mergers that form extreme 
helium stars or helium-rich sdO stars. 

Further progress requires advances in both theory and 
observation. Understanding how ELM WD binaries merge 
and ignite is an open problem. Greater numbers of AM 
CVn s tars must be fo und to better constrain their prop- 
erties (|Rau et alJl201Ch . New imaging surveys such as the 
Palomar Transient Factory, Pan-STARRS, Skymapper, and 
eventually the Large Synoptic Survey Telescope will discover 
large numbers of unusual transients such as underluminous 
SNe. 

Our goal is to use observations to investigate links be- 
tween double-degenerate mergers and unusual transients. In 
the future, obtaining a larger sample of 30 ELM WDs will 
improve the uncertainty in our ELM WD space density and 
merger rate estimates by a factor 2. We anticipate this is 
the first step in our ELM Survey towards a greater under- 
standing of merging ELM WD systems. 
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